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Abstract 

We examine a simple model of interaction of gravitational waves with matter (primarily 
represented by dust). The aim is to investigate a possible damping effect on the intensity of 
gravitational wave when passing through media. This might be important for gravitational 
wave astronomy when the sources are obscured by dust or molecular clouds. 

1 Introduction 

At the present time we witness an ongoing dedicated search for gravitational waves from various 
astrophysical sources using the ground based laser interferometers [1]. Since the sensitivity of 
these instruments is at the upper bound of assumed intensities coming from realistic sources it 
might be interesting to estimate the damping effects of various dust formations on the propagating 
gravitational waves. 

In the past several papers (for review see [2] ) dealt with the effects of media on gravitational 
wave but their main concern was the modification of dispersion relation and not the possible 
influence on the amplitude. Two different models for the medium most often used were: medium 
composed of deformable "molecules" with internal structure [2,3] giving rise to anisotropic pressures 
or free particles with rare collisions described by kinetic theory [4]. To describe the damping in 
the second model the rate of particle collisions has to be addressed giving rise to the imaginary 
part of the refractive index [5] . The damping was also studied in viscosity approximation [6] . In 
these treatments the damping was negligible due to the slow accumulation of the phase shift that 
is important for the damping effect only at the values close to 27r. 

In sections 2 and 3 the damping is derived using a straightforward computation without the 
assumption regarding particle collisions. The interaction of gravitational waves and matter (repre- 
sented by dust cloud) is derived based on the following picture. The incoming gravitational wave 
produces periodic oscillation within the "molecules" in the dust cloud computed using geodesic 
deviation equation (hence we use the first of the above mentioned models for the medium). In 
the second step these oscillations themselves produce gravitational waves derived using multi- 
pole approximation. Finally we compose them with the original waves propagating through the 
cloud from a distant source. We use the term dust because we assume no interaction between the 
"molecules" . It means that the model cannot describe effects connected with collective phenomena, 
on the other hand it can still accommodate selfgravitation of the dust cloud if it can be encoded 
in the background metric. 

In section 2 we assume that the wave is approximately planar at the region of the cloud and 
the background might be described using Minkowski metric so the source of gravitational wave 
has to be far enough. In section 3 we generalize the procedure to curved background and treat the 
example of a spherical shell of dust around central source of gravitational waves. To arrive at the 
solution in this case the generalization of quadrupole momentum tensor to this specific problem is 
given. Section 4 briefly summarizes results of a previous work [11, 12] and mentions its seemingly 



1 



unfavorable implication. However it is shown that the result concerns orders of magnitude that 
are not decisive for any concrete implications. 

2 Derivation for planar waves in Minkowski background 

We start with the simplest case where the geometry in the neighbourhood of the cloud is assumed 
to be flat. First, we recall the behavior of particle separations in the flat background (that we 
use for simplification in the region of the cloud) under the influence of gravitational wave (de- 
termined by geodesic deviation). However we will generalize the formula to the situation where 
the reference center of the "molecule" is geodesic while the particle (whose separation from the 
center we want to determine) is additionally sitting in a harmonic potential (described by "spring 
constant" A: > 0) centered at the initial position and feeling dissipative effect of its environment 
(described by "damping constant" 6 > 0). At the same time, this is the most straightforward 
generalization of the second order differential equation using only linear terms. Generalization of 
geodesic deviation equation to include nongravitational forces is usually called relative deviation. 
We suppose that the particle is in rest before the wave arrives and therefore only use the so called 
"steady state" solution of driven, damped oscillator equation that serves for generalization of the 
geodesic deviation equation 

4(r) = x^MS^ + iC/i-nr)]U C = , (1) 



where lo is the frequency of gravitational wave and i is a complex unit. The equation ([T]) expresses 
the position of particle B in proper reference frame of particle A lying at the origin. Transverse- 
traceless (TT) calibration was used to compute the Riemann tensor components needed in geodesic 
deviation equation. It should be noted that there exists TT system which agrees to first order in 
perturbation with the proper reference frame. In the following we will use exactly this system for 
further computations. 

Next, we need to know how these periodic movements generate gravitational waves. For this 
purpose we use the reduced quadrupole moment 




We compute this expression using the equation Dropping the particle labels and assuming 
that the density p of the dust stays approximately constant throughout the passage of the wave 
(cross-section changes only in the second order of perturbation) the result computed to linear terms 
in perturbation is the following 

^jfe = P / [(SjSfc - ^r^Sjk) + ^CixjX^K^-^ + Xkx^h^f - pmx'h-^f'Sjk)] d^x , (3) 
Jv 

where x denotes the positions before the wave arrival. It indicates, that generally, we cannot 
neglect the influence of matter on passing gravitational radiation, since the generated intensity 
depends on second time derivative of reduced quadrupole moment. We can see that this intensity 
would be produced by the terms in the second round bracket in equation ([3]). 
The generated perturbations gh in TT frame have the following form [9] 

ghf,it^r) = ^il{t-r) + 0{iy) , (4) 

where 

If, = P^IabPt - ^P^kP'^'^ab (5) 

and Pab = Sab — naUb is the projection operator to the subspace transverse to unit radial vector 
ria- 
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2.1 Example with cuboid 



Let us show the computation for one very simple situation. Suppose the cubical dust cloud has 
approximate sizes X^Y and Z (with respect to Cartesian coordinates) in the reference frame with 
origin in the clouds center and the z-axis pointing to the distant source of gravitational waves 
on one side and to distant observer on the other side. We assume that the wave from external 
source is approximately planar in the region of the cloud. We compute the second time derivative 
of highest order terms for the relevant components of reduced quadrupole moment in TT gauge 
(with — (dz)"') for a "molecule" (region where the driven, damped oscillator model described 
above might be applied) of linear dimension a obtaining 

il = iCpa^/i- (6) 

TT 

xy 



= 



where h'^{z = 0) ~ g-it^t jg planar perturbation in TT gauge from distant source (we have chosen 
the © polarization mode for simplicity). By superposing the generated gravitational waves ([4]) from 
all molecules in a given layer of dust cloud in the xy-plane we obtain a planar wave (neglecting 
the effects of finite dimensions of the cloud) 

where the planar density of "molecules" a can be computed as a^^ for a dense medium and equality 
= —ui'^h'^ was used. Now we add the generated wave (O to the original perturbation h"^ . This 
process will happen in subsequent layers with efficiency described by parameter k —< 0..1 > 
(roughly speaking it would happen in every 1/k slice) which gives the following exponential decay 
of the amplitude of gravitational perturbation as it travels through the cloud 

hA^hse-'^^^^^P''" , (8) 

where A is linear density of molecules in the direction of propagation and we have omitted coor- 
dinate indices, hs is the amplitude before the cloud and Ha after it. We may define a volume 
number density of "molecules" as p — aX. 

Since the exponential factor in ([5]) is complex the overall effect is twofold. The imaginary part 
is responsible for the phase shift and the real part for the possible damping. The real part of the 
factor can be written as — '^^'^^'■'^^ ujpZMa^ , where M is the mass of the "molecule". This number 
is non-positive and therefore really describes exponential decay of amplitude. 

Although this picture seems to be pretty naive there is direct correspondence between above 
mentioned second time derivatives of quadrupole moment and perturbations of stress energy tensor 
representing anisotropic stresses that generate gravitational waves via standard linearized Einstein 
equations (for detailed discussion see [2]). 



3 "Spherical waves" and the dust shell 

Now we are going to a more complex example where the approach from the previous section 
cannot be applied straightforwardly. As a background we use the Vaidya-(anti-)de Sitter metric 
in stereographic coordinates (u, r, 77, ^) 

= (_i + ^ H^r^)du^ + 2dudr + ^{dif + df) , (9) 

r p'^ 

where p{r],£,) = 1-1- 1/4(77^ -I- 1^^) and H = \/ A/3 (with A being a cosmological constant). To 
determine the form of gravitational waves we can use the results of paper [7] where the background 
is a general Robinson- Trautman spacctimc with cosmological constant (describing also the above 
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metric as a special case) and high-frequency approximation developed by Isaacson [8] is used. The 
wave vector is A:^ = {4>^m 0, 0, 0) and we specialize to the © polarization mode. So the only nonzero 
elements of the perturbation tensor are the following 

r 

Kv ^^J^ exp(i0(u)) = , (10) 

with A being an amplitude and </) the phase. The effect will be analyzed with respect to the 
stationary observer near infinity with four- velocity v = (zi, 0, 0, 0). For zero cosmological constant 
this observer is almost geodesic. Once again we will employ the geodesic deviation equation (in a 
more general form [9]) 

VvVvn" = -R'^.^^sv^n^v' . (11) 

The deviation vector rt" connecting neighboring geodesies has components in the dri and 5^ di- 
rections only. Assuming that the spherical dust shell that interacts with the gravitational wave 
starts the free fall from rest (^^ = 0) the right-hand side of pT|) could be written as '^^"2 and 
ii = (— ffuu)"^^^ from four-velocity normalization. The Riemann tensor could be decomposed into 
background part R^'^^ and contribution induced by perturbations i?*-^-* which is the dominant one 
in the high-frequency approximation [8] 

. d^ 

^^7^ - '^himo.m = -2~'^'^d^ exp(z</.(z.)) = ^ . (12) 

Now we can integrate the simplified equation (fTT|) (assuming that i?'-"-' changes very slowly) 
with added effect of non-gravitational forces (i.e., when the particle position satisfies VvVvH — 
— i?(v, n)v — foVvH — fc(n — no), generalizing the relative deviation used in the flat case) to obtain 

= 1 + + Z-^/^^ff"" exp(z0(u)) nl , (13) 



where C — 

Next we would need a fully relativistic generalization of the quadrupole moment defined using 
integrals over the source. However this seems to be a very tricky question and the only agreed 
definition relies on identifying the corresponding contribution in the multipole expansion of the 
field generated by the source which is not useful in our case. For example in the case of rigidly 
rotating axially symmetric star the multipole moments were numerically computed using iteration 
techniques and Green function method by Ryan [10]. But the method is too complicated for our 
purpose so we will try to generalize the flat-space formula since our two-space where the deviation 
vector resides is conformally flat (to remedy the problem with coordinates without simple length 
interpretation we use tensor with both types of indices). Our prescription for a small clement on 
the pole is the following {A, B take values 77 and ^) 



r2 



^A^P J [riAn'' - l/3r2,5f drdr^dC , (14) 

K{ro,ro+Ar)nAfl 

where Afl = ({— A77/2, Arj/2} x {— A^/2}), Ar] = = a/rg and the metric coefficient grr 
comes from Schwarzschild-like form of the metric. Integration is over the intersection of a shell 
between the indicated radii and small spatial angle. Using the integral mean value theorem, setting 
= gABn^ , expressing the deviation vector as Uq = rj {uq — ^) and using (|13p we obtain in the 
highest order (already using the TT-gauge) 

^ 2-3/2 + pACgrrArr'.-^ieM^c^iu))) J^^^ ^^d^d^ ^ (15) 



0-5/2 / JliO)v 

"""1 pACff„Arro-ia4(</),„)2exp(z<^(u)) -KO(a^), 
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where we have neglected the terms containing second derivative of phase in accordance with [8] 
and grr might be approximately set to one when we are not close to horizons. In the same way 

one obtains — '~Su^^'^v Sn^^'^v ~ ^' ■'^^'^i '^^ need to use equation ^ to obtain 

gravitational waves generated by the element and subsequently add the contributions from a given 
layer. In the spherical case we count (in integration) only those elements that are "visible" on 
the surface of the spherical layer from the point on the next layer (for the north pole it means 
9 ~< 0.. arccos( ^^^°^^ ) >), therefore not taking into account curved geometry where the waves 
propagate. When neglecting the term periodically dependent on the square root of distance (which 
corresponds to discarding the term related to boundary in the flat case) we obtain 



gii- n 



2mra ^ d^ ^tt») 
6 o,r dv? 



V 



(16) 



where a is the surface density of the "molecules" . For the covariant metric perturbation one may 
then write 

Kv (^0 + U + l)Ar) = [1 + grr\(r„+jAr}KAr] h^^{ro + jAr) , (17) 

where 

3 \ kguu 



K 



(18) 



— — — (with >> 1 being the number of layers and ro,ri the inner and 



N 



Setting Ar(=: A^^) 

outer radius of a dust shell, respectively) and assuming that the wave decreases with efficiency 
K(as in the planar model) we get for the intensities under inner radius hs and above outer radius 
hA 



N 
J=0 



KAr 



ro +jAr 



H^ro+jArf 



= hf 



1 K{N + 1)N {r^K - 2 m + 2 H^ro^) Ar^ 



(ro — 2 m 



{N + l)rQKAr 
ro — 2 m — H^vq^ 



O (Ar^) 



(19) 



We can check that this result has at least the same length dimensionality and similar dependence on 
frequency as for our cuboid case on Minkowski which might give some vindication to the calculation 
used in the curved space. 



4 Perturbation viewpoint 

In 1987 Ehlers, Prasanna and Breuer [11] investigated the propagation of small-amplitude gravita- 
tional waves through pressureless matter ('dust') using the linearization of Einstein equations for 
dust 

g^^, p, — > gf^„ + eSg^^, p + eSp, + eSUf_, . (20) 

They have used the WKB method to study the locally plane, linearized perturbations of an ar- 
bitrary background dust spacetime asymptotically for small wavelengths (thus assuming the high 
frequency character of the perturbations) 

%. = 5R[et^(^)/^,(a;,£)] . (21) 

Finally, they have cast the resulting system of equations into the form where Sg^i, serves as a 
source for generating perturbations Sp and Su^. 

However, they have arrived at the conclusion that in the leading order no perturbations of 
density or vorticity are generated 

Sp = Su^^O . (22) 
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This makes it impossible to alter the passing gravitational wave in the second step (the medium in 
this model is unchanged in the first place) . This result would not allow the damping of gravitational 
wave seen in the previous sections. But one should also keep in mind that the model from sections 
2 and 3 relies on an internal structure of the "molecules" (although their collection behaves globally 
like dust) and therefore is not directly comparable to the perturbation picture described above. 

Anyway, careful reading of the paper [11] (and of a more recent paper [12] dealing with the 
perfect fluid case) reveals that in fact the order in which the right-hand sides of equations ((22|) are 
supposed to vanish is But this is necessary condition for the perturbation method itself to be 
consistent and does not restrict existence of nonzero matter perturbations at the order e. 

5 Conclusion 

The above described simple model (presented in sections 2 and 3) might produce measurable effect 
on the intensity of detectable gravitational waves. The model assumes certain (simple) interaction 
of particles inside the "molecules" forming the cloud so it most closely resembles molecular cloud. 
The results described in section 4 were shown not to cause any important restriction on the model 
from the point of view of the perturbative approach. 

To give some notion of the strength of this effect we can compute the damping for some 
astrophysically relevant data. We use the flat model (section 2) with k = lu'^ (maximal damping 
effect, see ([T])), h = cj/lO (to stay safely at underdamped oscillations when the driving force is 
neglected), a = 10~^^m (typical small molecule dimension) and efficiency k — 10~^ (to account for 
limited internal structure). We will compute the length over which the amplitude of gravitational 
waves is damped by 10%. For dense medium (p ~ a^'^) this length is 10~^^pc for frequency of 
lOOHz (and 10""'^° for IHz). For average interstellar medium densities of our galaxy (IQ^m^^) the 
lenght is however 10^°pc for lOOHz (lO^^pc at IHz). 

The model is to simple and rough at this stage to make some final judgement considering the 
observational consequences. Some parameters are not set by clear physically founded estimates. 
The future work will concentrate on providing a more complex model capturing more precisely 
the structure of the medium. If one is concerned with molecular structures the approach should 
probably utilize quantum mechanics rather that classical one at the "microscopic" level. Also the 
treatment of the curved background should be more consistent (e.g. the model neglects the effects 
of curved background in the computation of perturbation produced by given layer). 
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